Quenched Slonczewski- Windmill in Spin- Torque Vortex-Oscillators 
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We present a combined analytical and numerical study on double-vortex spin-torque nano- 
oscillators and describe a mechanism that suppresses the windmill modes. The magnetization dy- 
namics is dominated by the gyrotropic precession of the vortex in one of the ferromagnetic layers. In 
the other layer the vortex gyration is strongly damped. The dominating layer for the magnetization 
dynamics is determined by the current polarity. Measurements on Fe/Ag/Fe nano-pillars support 
these findings. The results open up a new perspective for building high quality-factor spin-torque 
oscillators operating at selectable, well-separated frequency bands. 



The advent of spintronics lead to the development of 
exciting new concepts for nano-scale devices using the 
spin-degree of freedom of the electron besides its charge- 
property pp. Much effort has been spent on spin-torque 
nano-oscillators (STNOs) PHZ], which typically consist 
of two single domain ferromagnetic layers separated by 
a metallic spacer or a tunnel barrier, one with its mag- 
netization fixed (polarizing layer), the other one suscep- 
tible to torques (free layer). An electric current travers- 
ing the system perpendicular to the layers becomes spin- 
polarized and exerts torques on the magnetizations [5]- 
I10j . leading to magnetization dynamics of the free layer. 
These excitations are typically in the range of a few 
gigahertz and can be detected by measuring the time 
variation of the magnetoresistance (MR). The pinning of 
the polarizing layer can, for example, be achieved by ex- 
change coupling to an antiferromagnet ;3] or by extending 
its thickness and lateral dimension [2]. In the absence of 
pinning both ferromagnetic layers will be excited and in 
the case of increasingly symmetric STNOs, this results 
in a dynamic equilibrium state called the Slonczewski- 
windmill [8, 11 . In this state both layers' magnetizations 
rotate in the same direction with a constant relative an- 
gle, resulting in a vanishing MR time-dependence. 
Here we investigate STNOs containing two stacked mag- 
netic vortices, i.e., a system consisting of two ferromag- 
netic disks, each in a vortex state and separated by a 
metallic, nonmagnetic spacer. Employing analytical and 
numerical methods, we study the coupled spin torque- 
driven motion of the magnetizations in the two disks, 
which are not pinned by any of the above mentioned 
mechanisms. We find that in the double vortex system, 
Slonczewski-windmill modes are quenched by an intrigu- 
ing mechanism. Our results show that that the current 
polarity determines which disk is excited and thereby se- 
lects the STNO-frequency band. This property is shown 
to arise from a spin torque-mediated vortex-vortex inter- 



action. Thus, it is an entirely different principle than the 
spin accumulation based mechanism suggested by Tsoi 
et al. [12]. We analyze in detail the underlying torques 
and the resulting forces. The force exerted by one vor- 
tex onto the other can be split into two contributions; 
one part arising from the polarizer vortex in-plane mag- 
netization acting on the free vortex core (disk-core part) 
and another one which is due to the core-disk interaction. 
We compute the dependence of these terms on the lateral 
core-core distance. These results provide insight into the 
fascinating dynamics of coupled magnetic vortices. The 
theoretical findings are supported by our experimental 
data obtained from double- vortex Fe/Ag/Fe STNOs. 
The motion of the magnetic vortex in each of the disks 
is governed by the Thiele equation |13) which we write 
here for the vortex in the top disk, 
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G = —2iTK,(iiQM s L/j)e z is the gyro vector, where L 
and k are the disk thickness and the vortex core po- 
larity, respectively. X is the core position with re- 
spect to the disk center, W refers to the effective mag- 
netostatic potential in which the core is moving, and 
D = (anoM s /^)Lirln(R/ro) characterizes the damping 
of the vortex motion. The parameters R and Tq are the 
radii of the disk and the vortex core, and the indices 1 
and 2 correspond to the top and bottom disks, respec- 
tively. The spin-transfer torque-induced force acting on 
the vortex generated by a vortex-polarizer can be decom- 
posed into two contributions F x = Ff + F^. Ff arises 
from the in-plane magnetization of the polarizer and acts 
on the core of the free vortex. The second term Ff is 
caused by the core of the polarizer-vortex and acts on 
the in-plane magnetization of the free vortex. Both force 
contributions depend on the lateral core-core distance 
Following Ref. [T3], we obtain the expressions 
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FIG. 1. Spin-torque magnitudes exerted by disk 2 on the 
vortex in disk 1 and their dependence on the lateral core-core 
separation. The torque arising from the polarizer core, F\ is 
negligible at large distances, explaining the results obtained 
by micromagnetic simulations and presented in Ref. [15] . 
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where we introduce the vectors 
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Ci — ±1 defines the vorticity ("+" corresponds to coun- 
terclockwise, "-" to clockwise) and m 2 ,« refers to the ith 
component of the bottom vortex unit magnetization vec- 
tor. The top layer magnetization is written in spheri- 
cal coordinates (ip,9), where ip and 9 are the azimuthal 
and polar angles, respectively. In the calculation we as- 
sume rigid vortices with thickness-independent magneti- 
zation. The cylinder axis is chosen as the z-axis. Using 
the ansatz of Feldtkeller and Thomas [TB] for the out- 
of-plane core magnetization = exp (— a 2 r 2 ) with 

a 2 ss In 2/(25 nm 2 ) to mimic the experimentally obtained 
vortex core size in Fe [17] ) , we obtain the forces and their 
dependence on the lateral core-core distance I as shown 
in Fig. [T] For large I, the contribution of F^; to the 
total force becomes negligible. This explains the simu- 
lation results reported in Ref. [TS], where the influence 
of the polarizer-vortex core on the dynamics was found 
to be small. From Eq. (|3| we see that this is caused 
by the reduction of the function |Vy| with increasing 
distance from the top vortex core. In contrast to the 
asymptotic decrease of F^; to zero, the magnitude of Ff 



approaches a finite value for large I. For small distances 
I < 2tq w 10 nm, however, we observe that both torques 
fall to zero. This can be attributed to the gain in symme- 
try with decreasing core-core distance. The small torque 
introduced by F^ is neglected in our investigation on the 
dynamics of the system with the vortices coupled by the 
electric current. This is justified by the fact that we are 
interested in the general behavior of the solutions. The 
decrease of Ff at small / must however be included. The 
coupled Thiele equations read, with j := TijP/(4e) 
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The sign of j is positive for electron flow from the top 
to the bottom layer. For a quantitative analysis of the 
solutions, we use parabolic approximations to the effec- 
tive magnetostatic potentials. We let Wi(Xi)/|Gi| = 
6.28 ns~ 1 X 2 /2, resulting in a top vortex eigenfrequency 
of fi = 1.0 GHz. The bottom disk potential is fixed at 
W 2 /\G 2 \ = (5/3)WV|Gi| (/ 2 ° = 1.7 GHz). These fre- 
quencies are chosen to represent our Fe/Ag/Fe nanopil- 
lars with ferromagnetic layers with a thickness ratio of 
5/3. We use a current of j/|Gi| = (5/23) ns" 1 corre- 
sponding to about 1.12 x 10 12 A/m 2 (for P = 1), which 
is within the range of experiments. The spin-transfer 
torque-induced force F d is assumed to increase linearly 
from I — to I = 10 nm, from whereon it is set to the 
constant value F^ = 23 nm. 

The solutions are obtained numerically using Maple's 
rkf45 implementation |18j . The results can be summa- 
rized as follows: For positive currents and equal vortic- 
ities, the top vortex gyrates around the disk center on 
a trajectory of about 50 nm in radius, regardless of the 
core polarity. The sense of rotation is determined by 
the core polarization (counterclockwise for positive and 
clockwise for negative core polarity). The gyration fre- 
quency is 1.0 GHz. The bottom vortex adapts its fre- 
quency and sense of gyration according to the top vor- 
tex. A dynamic equilibrium develops with a constant 
phase difference between the vortices even in the case of 
opposite relative core orientation. This frequency adap- 
tion is accompanied by a strong reduction in radius of 
the bottom core trajectory: for parallel cores, the radius 
is about 0.7 nm while for the antiparallel configuration, 
the reduction is even more pronounced (approximately 
0.2 nm). For negative currents, the vortices switch roles: 
The bottom vortex gyrates on a large orbit (18 nm), while 
the top vortex trajectory is quenched (1.2 nm for paral- 
lel, 0.3 nm for antiparallel core alignment). In the dy- 
namic equilibrium the phase difference is constant and 
the gyration frequency is about 1.7 GHz. This frequency 
corresponds to the eigenfrequency of the bottom vortex 
while the sense of rotation is determined by its core po- 
larity. The gyration phase difference between the two 
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FIG. 2. Trajectory radius r and frequency / of the bottom 
vortex as functions of the phase <j> with respect to the top 
vortex that gyrates at the same frequency for positive (red) 
and negative (green) core polarity. The blue lines mark the 
eigenfrequency /° of the free running bottom vortex for pos- 
itive (b) and negative (c) core polarity. The shaded regions 
correspond to equivalent solutions but negative radius. 




5 10 15 20 25 30 



t[ns] 

FIG. 3. X-components of the top (red) and bottom (blue) 
core coordinates versus time for the case of symmetric disks. 
The cores are aligned parallel in (a) and antiparallel in (b). 
In the parallel windmill-mode only appears in a tran- 

sient time interval, but is hindered afterwards. The oscillation 
decays due to the low core-core separation and the related de- 
crease of the spin-transfer torque-induced force. 
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vortices depends on the relative core alignment. For pos- 
itive currents and parallel cores, the bottom core gyrates 
approximately 90° ahead of its top counterpart, while for 
antiparallel cores a 90° lag is observed. From Eqs. (5) 
and (6) it is clear that the solutions for opposite vortic- 
ities arc identical to those obtained for equal vorticities 
with a negative current polarity. 

For large enough \j\, the obtained characteristics of the 
dynamics are the generalization of the criterion found 
in Ref. 15 . In the model used by those authors, the 
polarizer was assumed to be a fixed, rigid vortex, and 
only magnetization dynamics in the other, free disk was 
allowed. In our case, both disks can be polarizing or 
free layer. For a given combination of vorticities C1C2 
and applied current polarity, the system responds with a 
damped and a dominant gyration, the former defining the 
polarizing and the latter the free disk. The current polar- 
ity determines which disk is dominantly excited. There- 
fore, the generalized jCC-criterion reads: For jC\C% > 
the top and for ]C\Ci < the bottom disk is excited. 
A general analytical solution for steady state trajectories 
is difficult to establish due to the complexity of the de- 
nominator in Eq. (4); however, it is instructive to analyze 
the situation when one trajectory radius is much larger 
than the other, a condition that, according to our numer- 
ical results, also holds in the vicinity of the limit cycle. 
We consider the case of positive j CiCg an d replace the 
denominator in Eq. (4) by d := \J Xf + Y± , yielding the 
following relations between the radius r of the bottom 
vortex trajectory, its phase <f> (relative to the top vortex) 
and the common frequency / of the two oscillators: 
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These relations are displayed in Fig. [2] and reproduce 
the behavior observed in the numerical solutions: For 
both cases of positive and negative bottom vortex core 
polarity, the bottom vortex can adapt to the (a priori 
arbitrary) frequency of the top vortex by adjusting 
the phase. Positive (negative) frequency corresponds 
to counterclockwise (clockwise) gyration. As displayed 
in Fig. [2ja), this phase shifting comes with a strong 
reduction of the orbit radius- or, in other words, a 
quenching of the windmill-modes. The dashed lines in 
the shaded regimes correspond to solutions of negative 
radius. Since a reverse of the sign of the radius is equiv- 
alent to a phase shift of n, these negative-r solutions 
are identical to the trajectories represented by the solid 
lines. By means of a phase adaption and reduction of the 
radius, the vortex can use a fraction of the spin-transfer 
torque-induced force to assist or counteract the force 
due to its magnetostatic potential. The resulting radial 
force component can differ strongly from the purely 
magnetostatic force. It may even lead to an inversion 
of the relation between the sense of gyration and the 
core polarity. A special case is the configuration, for 
which a frequency adaption is not necessary, i.e., if 
the two disks are identical and the cores parallel [Fig. 
3k)]. In this case, the vortices start rotating in phase, 
but as they reach the limit cycle and the core-core 
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FIG. 4. Resistance versus field for (a) I = —21 and (b) 
/ = +21 mA sample current. I < corresponds to electron 
flow from bottom to top electrode. Squares on the curves 
indicate the observation of magnetization dynamics for the 
corresponding field value, (c) Field dependence of the double- 
vortex mode frequency for the two electron flow directions as 
indicated. 

distance drops below 10 ran, the mutual spin-transfer 
torque-induced force decreases leading to a decay of 
the oscillation amplitudes. For antiparallel alignment, 
the windmill-modes are quenched by the mechanism of 
frequency and phase adaption [Fig. [3Jd]. 
For an experimental confirmation of the frequency and 
phase adaption mechanism and the related quenching 
of the windmill modes we study the current-induced 
magnetization dynamics of an Fe/Ag/Fe nanopillar with 
a Fe layer thickness ratio of 5/3. According to our model, 
we expect to observe excitations for both current polari- 
ties, but with different frequencies yielding a frequency 
ratio of approximately 5/3. Cylindrical nanopillars 
are patterned using e-beam lithography and Ar ion 
milling from molecular beam epitaxy-grown GaAs(OOl)/ 
Fe(l)/Ag(150)/Fe(25)/Ag(6)/Fe(15)/Au(25) stacks 
(layer thicknesses in nm). The pillar diameter is 210 nm. 
The milling was stopped after reaching the 150 nm thick 
Ag buffer layer. Thus, the oscillator consists of two 
ferromagnetic disks of equal diameter and comparable 
thickness stacked on top of each other. Figures |4ja) 
and (b) display the field dependence of the nanopillar 
resistance for / = +21 mA (+6.1 x 10 7 A/cm 2 ), re- 
spectively. The external magnetic field was applied in 
the sample plane. The magnetoresistance profiles are 
characteristic for this sample type [TH] and reflect two 
magnetization states: The first one comprises a vortex 
in one disk, while the other nanomagnet remains in 
a quasi-homogeneous state. These configurations are 
characterized by a nearly linear field dependence of the 
resistance caused by a continuous lateral displacement 
of the vortex with changing field. The second state is 
observed from low field magnitudes up to about 200 mT 



and is characterized by low resistance values near the 
level in magnetic saturation. Here, each disk contains 
a vortex with the vorticity given by the circumferential 
Oersted field. This results in locally parallel alignment 
of the two disks' magnetizations explaining the observed 
low resistance. For both current polarities in Figs.[4ja) 
and (b) we detected magnetization dynamics in those 
field intervals, in which the double-vortex state occurs. 
The excitation frequencies are shown in Fig. |4jc) 
along with the corresponding electron flow directions of 
the externally applied currents I. All frequencies are 
below 2 GHz, which is typical for vortex gyration in 
Fe/Ag/Fe nanopillars [19 - 21 , but the frequencies are 
clearly different and well separated for the two current 
polarities. At low external fields their ratio is about 1.46. 
Using previous data on single vortex dynamics [THl [22] , 
we estimate the influence of the Oersted field on the 
zero-field vortex frequencies to be ~ 145 MHz for both 
disks. This shifts the above value to 1.54, i.e. very close 
to the ratio of the disk aspect ratios (5/3 ~ 1.67). Since 
the gyration frequency of a vortex in a ferromagnetic 
disk is in a first approximation proportional to the 
disk aspect ratio [33], our data strongly suggests that, 
as proposed, the observed signals at opposite current 
polarities originate from excitations in the two different 
disks, the low(high)-frequency mode corresponding to 
gyrotropic motion of the vortex in the top (bottom) 
disk. The dominantly excited disk is determined by 
current polarity in agreement with our model. 
In summary, we have theoretically and experimentally 
demonstrated that windmill modes are quenched in 
double-vortex spin-torque nano-oscillators. The origin is 
frequency and phase adaption of the gyrotropic motions 
in the two disks, which results in a strong suppression 
of the gyration radius in one of the disks. Changing the 
sign of the exciting current provides an effective mode 
selection mechanism, which allows to deliberately choose 
between separated frequency bands of the oscillator. 
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